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ABSTRACT 



in : 

i Aims. An all-sky survey of loop- and arc-like intensity enhancements has been performed in order to investigate the large-scale 
, structure of the diffuse far-infrared emission. 
_^ Methods. We used maps made of 60 and 100 fim processed IRAS data (Sky Survey Atlas and dust infrared emission maps) to 
\Q , identify large-scale structures: loops, arcs or cavities, in the far-infrared emission in the Galaxy. Distances were attributed to a 
f*"^ • subsample of loops using associated objects. 

Results. We identified 462 far-infrared loops, analyzed their individual FIR properties and their distribution. This data forms 
the Catalogue of Far-Infrared Loops in the Galaxy. We obtained observational estimates of /i n 30% and / out » 5% for the hot 
1 gas volume filling factor of the inward and outward Galactic neighbourhood of the Solar System. We obtained a slope of the 
O power law size luminosity function (3 = 1.37 for low Galactic latitudes in the outer Milky Way. 

Conclusions. Deviations in the celestial distribution of far-infrared loops clearly indicate, that violent events frequently overwrite 



Or 
i 

o 

^2 the structure of the interstellar matter in the inner Galaxy. Our objects trace out the spiral arm structure of the Galaxy in the 
. . neighbourhood of the Sun and their distribution clearly suggests that there is an efficient process that can generate loop-like 
>>• features at high Galactic latitudes. Power law indices of size luminosity distributions suggest, that the structure of the ISM is 
ruled by supernovae and stellar winds at low Galactic latitudes while it is governed by supersonic turbulence above the Galactic 
plane. 
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1. Introduction 

The large scale structure of the cold interstellar matter 
can be significantly affected by violent events. This struc- 
ture - which is diverse with the complex distribution of 
shells, cavities, filaments, arcs and loops - is often re- 
ferred to as the "Cosmic Bubble Bath" (Brand & Zealey 
I1975[) . The evolution of bubbles and superbubbles pro- 
duced by supernova explosions and stellar winds of as- 
sociations are the primary processes that determine the 
structure and energetics of all components of the diffuse 
interstellar medium. Studying supernova (SN) explosions 
in a uniform medium, Cox and Smith (1974) pointed out 



that if the Galactic SN rate is "sufficient" it can produce 
the " swiss cheese morphology" of the cold diffuse ISM with 
hot coronal gas inside the bubbles. This was built into the 



model by McKee and Ostriker (1977) where SN explo- 
sions in the cloudy ISM produce a three-phase medium. 
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In this scenario the next generation of stars are born in 
the compressed medium of the bubble walls (propagating 
star formation, see Blaauw 1991, for a review). 

The idea that supersonic turbulence may be impor- 
tant in the control of star formation became recognized in 
recent years (see MacLow & Klessen [20041 for a review). 
Turbulence - if dominant can also rule the structure of 
the ISM creating clouds and cavities with regular proper- 
ties of a fractal geometry, apparently similar to the struc- 
ture formed via violent events. 

Structures produced by various processes are often 
characterized by the power law parameter f3 of the me- 
chanical luminosity distribution of bubbles or holes. With 
a specific model this can be converted into a parameter 
describing the size distribution of these structures, then 
[3 can be determined even for a medium with no stellar 
energy injection (e.g. a str ucture created by supersonic 
turbulence). Oey & Clarke |( 1997)1 obtained 1.4</3<2.1 
for the nearby galaxies M31, M33, Holmbergll, SMC. 
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Ehlerova and Palous (2005) deduced /3=1.6±0.3 from the 
size distribution of HI shells in the outer Milky Way. 
Elmegreen (|1999j) derived /? = 2.15 for a medium with a 
structure of D = 2.3 fractal dimension, also from the size 
distribution of holes. This structure should be representa- 
tive of that created by supersonic turbulence. 

The hot gas volume filling factor / of a galaxy is 
also an important parameter in characterizing the lifecy- 
cle of the interstellar medium. Although there are several 
estimates of / based on theoretical considerations and 
semi-empirical models, the / parameter of our Galaxy 
is hardly constrained observationally (Ferriere 1998 and 
Gazol-Patino & Passot 1999). Recently Ehlerova and 
Palous ( 2005) obtained /=5% for the shell volume filling 
factor for the outer Galaxy through an automated identi- 
fication of HI shells. 

To derive these global parameters a large sample of 
object are needed, and it is not easy to construct such a 



database for our own Galaxy. Thilker et al. (1998) and 
Mashchenko et al. (1999) made efforts to detect HI shells 
automatically, based on a model of the shell. Recently 
Ehlerova et al. ( 2004) developed a model-independent al- 
gorithm which automatically searches HI shells in data 
cubes and identified ~1000 structures in the Leiden- 
Dwingeloo Survey data. An artificial neural networks al- 
gorithm by Daigle et al. (2003) was successfully applied to 
the Canadian Galactic Plane Survey data. However, these 
works intended to detect HI holes based on the velocity 
information of the 21cm line data rather than on mor- 
phology, therefore they cannot find bubbles at late stages 
of their evolution when the expansion has slowed down. 

Shell- or arc-like intensity enhancements are reported 
in many tracers of the ISM (see Kiss et al. 120041 and ref- 
erences therein, for an introduction) . The structure of the 
diffuse ISM is well-represented in the far-infrared, mainly 
observed as the Galactic cirrus emission (Low et al. 1984). 
Although some studies reported loop features identified 
on the far-infrared, these were either restricted to the 
Galactic midplane (Schwartz 1987), performed for a spe- 
cial object type, e.g. Wolf-Rayet stars (Marston |l996). 
Most of the prominent HI loops are conspicuous in the far- 
infrared, too, due to their dust content, which is closely re- 
lated to HI (Boulanger & Perault |T988) . Many individual 
loops of various origin were first identified on FIR images. 
E.g. the North Celestial Pole Loop ( Meverdierks 1199 f|) is 
an example of a high velocity cloud Galactic disc in- 
teraction. The Cepheus Bubble is a remarkable infrared 
loop l)Kun et al. 1987(1 which is likely a combined result 
of multiple supernova explosions and stellar winds in the 
Cep OB 2 association. 

In a recent work Kiss et al. (2004, hereafter KMT04) 
presented the results of a quest for far-infrared loop fea- 
tures in the 2 nd Galactic Quadrant. They catalogued 145 
loops and investigated their morphological and physi- 
cal characteristics. This study was not restricted to the 
Galactic midplane, but reached high Galactic latitudes 
and was able to give a comprehensive view on the dis- 
tribution of large scale intensity enhancements. 



In this present work - as a continuation of the 
KMT04 study - we extend the quest to the 1 st , 3 rd and 
4 th Galactic Quadrants. These two works together provide 
the Catalogue of Far-InfraRed Loops in the Galaxy. 

2. Input data and data analysis 

We used the 60 and 100 ^m IS SA plates (IRAS Sky Survey 
Atlas, Wheelock et al. 1994) in order to explore the dis- 
tribution of dust emission. Reprocessed IRAS 100 /im and 
reddening maps by Schlegel et al. (1998) were investigated 
as well. We searched for loop- or arc like features in these 
images in the 1° < D < 40° diameter range. The data re- 
duction steps and derived parameters are the same as de- 
scribed in KMT04. For the detailed description of the data 
reduction and derived parameters we refer to KMT04 and 
Appendix A (online only) of this present paper. 

3. Results 

In our all-sky survey 462 FIR loops have been iden- 
tified. We found 317 loops in the 1 st , 3 rd and 
4 th Galactic Quadrants (presented in Appendix B, on- 
line only), beyond the 145 loops in KMT04 (2 nd Galactic 
Quadrant). The electronic version of the catalogue 
with additional data products can be found at: 
" http: / /kisag.konkoly.hu/CFIRLG" . 

The catalogue contains the name of the loop; the cental 
Galactic coordinates, size and position angle of the fitted 
ellipse; and other parameters describing the appearance 
of the loop in FIR images. For details, see Appendix B in 
the electronic version of the paper or Sect. 3 in KMT04. 
For a subsample of 43 loops of our catalogue (1 st , 3 rd and 
4 th Galactic Quadrants) we were able to derive distances, 
using the distances of associated objects. The selection 
criteria for distance indicators were the following: 

i) objects related to possible energy injection sources 
(O/B stars, SN-remnant, etc.) had to be placed in the 
interior of the loop; 

ii) density enhancements of the ISM (molecular /dark 
clouds, molecular cores, etc.) had to be in the interior 
or in the loop wall; 

iii) there had to be at least two objects with similar dis- 
tances; 

iv) distances of individual objects had to agree within the 
uncertainties. 

In many cases the distance of a loop has already been de- 
termined by previous studies. In Appendix C we present 
the distances and a detailed list of distance estimators for 
these 43 loops. Fig. shows the distribution of estimated 
distances and the distances projected to the Galactic 
plane. 

In most cases, however, no definite distance could be 
derived for a specific loop. In these cases we listed all pos- 
sible associated objects (not taking into account the dis- 
tances of the individual objects), which appear on the loop 
wall or in the interior. The list of these associated objects 
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Fig. 1. Distribution of GIRLs in the sky (Galactic coor- 
dinate system, Mercator projection), represented by the 
fitted ellipses. Note that Mercator projection causes a size 
distortion at polar regions. 

is presented in Appendix D (online only, available at CDS: 
httpTA-.). 



4. Discussion 

4.1. Sky distribution 

The Galactic longitude distribution of the GIRLs partly 
reflect the spiral structure of the Galaxy (Figs. |3 and 
|3J). There is a noticeable increase in loop counts towards 
the major spiral arms i.e. the Local-, the inner Carina- 
Sagittarius- and the outer Perseus Arms. 

The distribution of GIRLs in the sky is expected to 
reflect the exponential disc distribution of the ISM, ac- 
cordingly regions closed to the Galactic plane should be 
the most populated parts of sky in loops and this num- 
ber should decrease nearly exponentially by the increasing 
Galactic latitude. 

The distribution of loop centres in Galactic latitude 
(Fig.QJ - which is expected to reflect the exponential disc 
distribution of the ISM - shows relatively large counts at 
high \b\ values, especially for outward locations. This can- 
not be explained by projection and distance effects only, 
in agreement with that was found for the 2 nd Galactic 
Quadrant (KMT04). To explain the formation of loops 
at high Galactic latitudes one needs an efficient process 
which does not necessarily have to take place close to the 
Galactic plane. This excludes SN-explosions and the stel- 
lar wind of massive stars as dominant effects. Although 
clouds infalling from the Galactic halo could create loops 
above the midplane, their infall rate is insufficiently small 
qFhlerovafo Pa Ions 199611. 

4.2. Groups at larger scales 

The celestial distribution of the identified loops is rather 
complex showing structures with scales even larger than 
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Fig. 2. Distribution of loops in Galactic longitude. The 
peaks in loop-count reflect the higher number density 
in the direction of Galactic arms (OS: Orion Spur, PA: 
Perseus Arm, C-S: Carina-Sagittarius Arm). The binsize 
is 10°. 



the average loop diameters (Fig. [TJ. E.g. two remarkable 
features, forming 'chains' of loops can be found around 
Z«175°, 6«+10° and around I « 115°, &w-15 Q . This lat- 
ter huge structure is well seen as a fragmented arc in HI 
too, at velocities approximately -40 kms^ 1 in the Leiden- 
Dwingeloo Survey (Hartmann & Burton 1997). 

Some groups of loop centres coincide well with the lo- 
cation of molecular complexes, including the Ophiucus 
Molecular Cloud (Z«10°, bw+15°), the Ursa Maior 
Molecular Cloud (Z«135°, 6«+55°), the Taurus molecu- 
lar complex (i«180°, 6w+0°) and the Orion-Monoceros 
molecular region (Jw210°, bw-15°). 

Although it is likely, that loop structures are associ- 
ated with large scale molecular material, no other obvious 
occurrences were found on these scales in our sample. 

4.3. Dominated by confusion? 

To independently study the immediate inner and outer 
surrounding of the Solar System, we arbitrarily split the 
Galactic latitude space into two parts, approximately 
along the direction of the Local Arm, i.e. the / = 70°-250° 
axis. The loop counts in the inner and outer part show sig- 
nificant differences: there are almost twice as many loops 
in the outer region than in the inner one. One possibility 
to explain this fact is the existence of strong confusion in 
the direction of central (inward) Galactic regions, as dis- 
cussed below (for the definition of Galactic regions A, B, 
C, D see Table 1). 

The possibility of a confusion-effected distribution for 
the 2 nd Galactic Quadrant was already mentioned in 
KMT04. In the present study - due to the extension to 
the whole sky - this confusion effect is more expressed and 
manifests itself in two main ways : 
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Region 


Galactic Longitude 


Galactic Latitude 


A 


I < 70° or I > 250° 


\b\ <30° 




inner 


low 


B 


70° < / < 250° 


|6| <30° 




outer 


low 


C 


I < 70° or I > 250° 


30°<|6|<60° 




inner 


high 


D 


70° < I < 250° 


30°<|fe|<60° 




outer 


high 



Table 1. Definition of four Galactic regions. Regions have 
been defined to help the comparison of loop characteristics 
at different parts of the sky. 



i) As presented in Fig. we see more loops in 
the 'outer' Galaxy than in the 'inner' regions for 
|6| < 30° (region A), although probabilities for the fre- 
quency of high-pressure events (e.g. SN-rates) should 
be significantly higher in the inner Galaxy. 

ii) The size distribution is different for regions A and B 
(see Fig. 03): in the inner Galaxy (A) the average appar- 
ent size is smaller than in B (4?1 and 5? 2, respectively). 
This discrepancy is not observed for higher Galactic 
latitudes (C and D). 

These two facts together indicate that (a) either we miss 
large loops at low \b\ in the inner Galaxy due to the strong 
IR background (small loops are still visible since they are 
smaller than the characteristic scale of molecular clouds) 
or (b) large loops are really destroyed by frequent vio- 
lent events and are only visible at their early evolutionary 
phases (i.e. at small size). Unfortunately our present FIR 
data cannot distinguish between these possibilities. 

4.4. The large scale distribution of ISM in the Galaxy 

In Appendix C we presented a subsample of 43 loops of our 
catalogue. For these loops we were able to derive distances 
using the distances of associated objects. Together with 
the 30-loop sample in appendix B in KMT04 we have a 
sample of 73 loops with known distances. We used this 
sample to characterize the large-scale distribution of the 
ISM in the vicinity of the Solar System. 

4.4.1. Size distribution of cavities 

In Fig. we present the differences in distance distribu- 
tions in the inward and outward Galactic regions. In the 
inward region A most of the loops are squeezed into a 
small area around the Solar System, within a distance of 
~0.2kpc (see Fig. 7). The rest of the loops do not show 
a concentration at a specific distance but are rather ran- 
domly distributed in the interval 0.4 kpc < d < 4 kpc. Note, 
that the 2 nd Galactic Quadrant has a major contribution 
to outward regions, therefore this distribution is quite sim- 
ilar to the one found in KMT04. 

We checked the distances of the individual loops (those 
with known distances in regions A and B) and the dis- 



tance of the wall of the Local Bubble in the direction of a 
specific loop using the maps and data based on the mea- 
surements of the Nal D-line doublet by ILallement et al.l 
(2003). The reassuring result of this comparison was, that 
all of these loops were in or behind the wall, i.e. behind the 
W(D2) = 20 mA contour of equivalent width (see Fig. |HJ. 

The general structure of galaxies is often character- 
ized by the distribution of bubbles and super bubbles , as 
was done for nearby galaxies by Oey & Clarke (1997) and 
Kim et al. |2003). The generally used parameter /3, the 
power-law index of the mechanical luminosity function, is 
closely related to the power law index s of the logarithmic 
size distribution: s — 2-2/3 (see the model by Oey & Clarke 



(1997) I. Using this relationship we derived a power law in- 
dex of (3 = 1.37 ± 0.17 for our outer sample (region B) for 
the interval 1.4 < log(2i? e ff) < 2.6, which is similar to the 
value of the 2 nd Galactic Quadrant only (f3 = 1.24 ±0.30, 
KMT04) and is in good agreement with (3 — 1.6 ± 0.3 ob- 
tained by Ehlerova & Palous (2005). Our value is lower 
than that of the other investigated galaxies (see Fig. ^jp) . 
The double-peaked distribution observed in KMT04 is less 
expressed in this larger sample and is closer to a power- 
law. Due to the insufficient sample no proper j3 could be 
derived for the inner Galaxy (region A), however, the dis- 
tribution is very shallow (Fig. 9a). 

We have no information on the distances of loops at 
high Galactic latitudes due to the lack of proper distance 
indicators. This prevented us from deriving proper (3 val- 
ues for high Galactic latitudes (C and D) as we did for A 
and B. However, the distribution of apparent size of the 
known-distance sample and that of the full high Galactic 
latitude sample are remarkably different. This indicates 
different real size distributions behind the apparent ones, 
and the presence of a different process governing the struc- 
ture. 

Due to the lack of distance information only the ap- 
parent size distribution can be used as diagnostic tool to 
test the structure of the high Galactic latitude ISM. In 
a simple model we tested the apparent size distribution 
with different [3 values. The model included a real size 
distribution with a prescribed (3 value (so that the slope 
of the log(ii) vs. log(iV[log(i?)]) distribution is s = 2-2(3) 
and a certain distribution of cavity centres. We derived 
the apparent size distribution from this size/spatial loca- 
tion distribution. The apparent size distribution was then 
fitted by a power law (as in the case of the real size dis- 
tribution), with a power law index of s Qpp , in the interval 
1° < 2i? app < 40°, i.e. in our investigated apparent diame- 
ter (2i? app ) range. 

The main features and results of this simple model 
were the following: 

— Assuming a random, homogeneous distribution of cav- 
ity centres, with a minimal distance of d m j n = pc, the 



final 



>app 



has little dependence on the absolute value 



of the maximal distance d max , as long as the appar- 
ent sizes are in the fitted range. Realistic estimates of 
the maximal distances are from ^250 pc (scale height 
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of the Galactic disc; Nakanishi & Sofue, 2003) to a 
few kpc (the maximal visibility distance), depending 
on the Galactic region we try to model. 

— The assumption of d m i n = is not realistic, since the 
Solar System is located inside the Local Bubble, and 
the loops must be behind its wall, as it is the case for 
the objects in our known-distance sample. We choose a 
minimal distance of d m i n = 130 pc (a characteristic dis- 
tance of the wall of the Local Bubble, see Lallemcnt 
et al. 2003), while keeping the random, homogeneous 
distribution behind. In this configuration the appar- 
ent size distribution is clearly altered, compared to 
the d m i n = cases. The apparent size distribution is 
shallower, and the larger the maximal distance, the 
closer the s app values with d m j n ^ to the ones with 
d m in = 0. Although the difference is noticeable, s app 
has a stronger dependence on (3 then on d m ;„ (not con- 
sidering very extreme conditions). 

— We have checked with a Kolmogorov-Smirnov (KS) 
test, whether the apparent size distribution of high la- 
titude loops in the northern and southern hemispheres 
may originate from the same distribution. The KS 
test resulted in a significance level of ~50%, show- 
ing that the two distributions are not really different. 
The Sapp. high values derived independently for the two 
hemispheres resulted in very similar values, therefore 
we used the data of the two hemispheres together in 
the following. There is a noticeable difference in the 
total counts (110 and 78 at the northern and southern 
hemispheres, respectively), but this does not affect the 
Sapp, high values. 

— Assuming /3=1.37 (that of the real size distri- 
bution in the known-distance sample), the ap- 
parent size distribution has s app =-0.7 for the 
d m in = 0, and -0.7 < s app < -0.5 for the d m i n = 130pc 
case, choosing the maximal distance in the range 
250 pc < d max < 1 kpc. This shows a good agreement 
with the power law index of the apparent size dis- 
tribution of the known-distance sample, s a p Pj known — ~ 
0.67±0.12. 

— The known- distance sample is characteristic for low 
Galactic latitudes (b<30°) only. The power law in- 
dex of s app , high =— 1.78±0. 21 is quite different from 
s app , known- The likely difference is a different real size 
distribution, i.e. a different (3. We tested, whether a re- 
alistic (3 can preproduce this value: /3 = 2.15 (s = -2.3, 
the canonical value for a fractal structure created by 
supersonic turbulence) resulted in s app =-1.8 for the 
d m in = and -2.2 < s app < -1.8 for the d min = 130pc 
case. These values are indeed very similar to s apPi hi g h- 

— We checked whether a size distribution with (3 = 1.37 
may be responsible for an apparent size distribution of 
Sapp = —1-78, assuming a specific spatial distribution 
of cavity centres. To reproduce the required values of 
both parameters at the same time, the spatial distri- 
bution has to be so, that the loops are highly concen- 
trated at the outer edge of the investigated volume and 
practically no loops are found at small distances. This 



is an unlikely configuration even taking to account that 
the Solar System is located inside the Local Bubble. 

As we saw, various (3 values lead to significantly differ- 
ent apparent size distributions, while changing the mini- 
mal and maximal distances in a realistic interval alters the 
size distribution less. Therefore it is not very likely, that 
very different (3 values would lead to the same s app us- 
ing different spatial distributions, in the framework of our 
model. We have to note, that among other limitation like 
the unknown real spatial distribution, our simple model 
do not take into account any observational biases, which 
are especially important at small apparent sizes. 

However, our simple model clearly outlines the differ- 
ence in the real size distributions at low and high Galactic 
latitudes. A probable explanation for this difference is that 
the structure is governed by different physical processes. 
Our high latitude apparent size distribution favours a (3 
value similar to that of supersonic turbulence (Elmegreen, 
1999), which has been long considered as the most likely 
process forming the structure of the interstellar medium 
further from the Galactic midplane. 

4.4.2. Volume filling factor 

The hot gas volume filling factor / is predicted to change 
with the Galactocentric distance (Ferriere|T22HIand Gazol- 
Patino & PassotCnnni: it is 20. ..30% for the inner Galaxy 
and decreases significantly for larger Galactocentric dis- 
tances. For low Galactic latitudes (&<30°) - where the 
structure is mostly formed by violent events - one can as- 
sume, that the distribution of hot gas is very near to the 
distribution of holes/cavities and therefore can be charac- 
terized by their distribution. Applying the same procedure 
as in KMT04 (see Sect 4 in KMT04, paragraph "Volume 
filling factor") for the whole loop sample in the inward 
and outward regions separately for b < 30° (regions A 
and B), we obtained /i n ~30% and / out ~5% for A and 
B, respectively. This latter value is very similar to the 
one presented in KMT04 for the 2 nd Galactic Quadrant 
(4.6% < /kmto4 = 6.4%). The size distribution in the in- 
ward region A almost completely misses the large loops 
which might be due to confusion. However, these may not 
exist at all, since large structures are regularly destroyed 
by violent events. The /; n value we found is in the or- 
der of ^20%, as predicted by Fcrricrc (1998) and Gazol- 
Patino & Passot l|1999f) . Our / out value is very similar to 
/hi ~ 5% found by Ehlerova and Palous (2005) through an 
automated identification of HI shells in the 2 nd Galactic 
Quadrant. 

5. Summary 

We performed a survey of loop/arc-like intensity enhance- 
ments in the diffuse far-infrared emission in the Galactic 
longitude intervals 0° < I < 90° and 180°</<360°. 
Merged with the results of the second Galactic Quadrant 
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(KMT04), we altogether identified and catalogued 462 of 
these features. 

The Catalogue of Far-Infrared Loops in the Galaxy con- 
tains the basic physical properties and a list of associated 
objects as well. The electronic version of the catalogue is 
available at: http://kisag.konkoly.hu/CFIRLG/ We also 
gave distance estimates for 73 loops based on the distance 
of associated objects. Our database provides a great op- 
portunity to study the large scale structure of the ISM in 
the Galactic neighbourhood of the Sun. We determined 
observational estimates for the hot gas volume filling fac- 
tor of the the inner (/j n ) and outer (/out) Galactic envi- 
ronment of the Solar System and obtained f- m ~ 30% and 
/out ~ 5%. These values are in good agreement with both 
theoretical models and values derived from measurements 
in the radio domain. A study of the cavity size distribution 
has shown that the structure of the ISM is different at low 
and high Galactic latitudes. While at low latitudes it can 
be explained by high pressure events, the apparent size 
distribution of holes favours a fractal structure at high 
latitudes, similar to what can be created by supersonic 
turbulence. 
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Appendix A: Detailed description of data 
reduction steps and main parameters 

Identification of the loop structure: We inves- 
tigated the 60 and 100 ISSA plates (IRAS Sky 
Survey Atlas, Wheelock et al. 1994) to explore the dis- 
tribution of dust emission. We created composite im- 
ages of the 12?5 x 12?5 sized individual ISSA plates us- 
ing the " geom" and " mosaic" procedures of the IPAC- 
Skyview package 1 , both at 60 and 100 /an. These im- 
ages were built up typically from 10. ..15 ISSA plates, 
reaching a size of ~ 40° x 40° . Loop- like intensity en- 
hancements were searched for by eye on the 100 ^m 
mosaic maps. Loops by our definition must show an 
excess FIR intensity confined to an arc-like feature, 
at least 60% of a complete ellipse-shaped ring. A loop 
may consist of a set of bright, more or less isolated, 
extended spots, or may be a diffuse ring or part of a 
ring. The size of the mosaic image limits the maxi- 
mal size of the objects found. On the other hand, due 
to the relatively large size of the investigated regions, 
loop-like intensity enhancements with a size of <1° 
were not searched for. The original ISSA /|o SA and 
/Jqo a surface brightness values were transformed to the 
COBE/DIRBE photometric system, using the conver- 
sion coefficients provided by Wheelock et al. (1994): 

- I 60 = 0.87 x 7|SSA + o 13 MJysr -i 

- Jioo = 0.72 x /issA _ L47 M j ysr -i. 

Dust IR emission maps by Schlegel et al. (1998) (SFD) 
were investigated to derive parameters describing our 
loop features (see Sect. 2.2). The main differences of 
the SFD 100 /zm map with the ISSA maps are the fol- 
lowing: (1) Fourier-destriping was applied; (2) aster- 
oids and non-Gaussian noise were removed; (3) IRAS 
and DIRBE 100 /im maps were combined, preserving 
the DIRBE zero point and calibration; (4) stars and 
galaxies were removed. 

We analysed the radial surface brightness profiles of 
the loops on the SFD 100 ^m map in order to check the 
effect of the removal of the sources mentioned above. 
We also used the SFD E(B-V) maps derived from the 
dust column density maps. In the case of these maps 
the colour temperature was derived from the DIRBE 
100 and 240 /im maps, and a temperature-corrected 
map was used to convert the 100 ^m cirrus map to a 
map proportional to dust column density (see SFD). 

— Ellipse fitting: The shape of a loop candidate was 
approximated by an ellipse, which was then fitted us- 
ing a 2D least-squares fit method. The fitted ellipse is 
defined with the central (Galactic) coordinates, the mi- 
nor and major semi-axis of the ellipse, and the position 
angle of the major axis to the circle of Galactic latitude 
at the centre of the ellipse. This latter was defined to 
be "+" from East to North (or counter-clockwise). 

— Intensity profile: For all of our loops we extracted 
radially averaged surface brightness profiles, extending 



1 http:/ /www.ipac. caltech.edu/Skyview/ 
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to a distance of twice the major (and minor) axis of 
the fitted ellipse, using 40 concentric ellipsoidal rings. 
These surface-brightness profiles (ISSA 100 and 60 fj,m, 
SFD 100 /jm and SFD reddening maps) were used in 
the following to determine the basic parameters of the 
FIR emission in the loop; two examples are presented 
in Fig. A.l. 

— Parameters of the loop wall: A local background 
was determined using the "non-loop" points in the 
radial surface brightness profiles, fitting a 3rd order 
polynomial, as demonstrated in Fig. A.l. This appro- 
priate background was then removed from each sur- 
face brightness profile point. The intensity excess I ex 
was derived as the maximum value of this background- 
removed profile. Inner and outer edges of the loop 
wall along the major axis a m and a out , respectively, 
are defined as the radial distance at the full width at 
/ex/2, the half power of the background-removed in- 
tensity profile. We distinguish three regions for a spe- 
cific loop: (1) loop interior (a < ai n ); (2) loop wall 
(a; n < a < a out ); (3) outer region ( a > a out ). The 
relative width of the wall of the fitted ellipse is defined 
as W = 1 - ain/oout- 

— Significance: We calculated the standard deviation 
of the background-removed surface brightness, c ex , in 
the 'non-loop' positions (a<ai„ and &>& aout ), and 
defined the significance of the loop as * = I cx /cr cx . We 
derived significance parameters on 100 and 60 /im ISSA 
maps and on the SFD 100 /im point source removed 
sky brightness and reddening maps (^noo, ^160) ^sioo 
and vI/sebv, respectively). The higher the value of 4" 
the higher the intensity excess of the loop over the 
background, i.e. ^ can be used as a quality indicator. 
It was a requirement, that a loop in the catalogue must 
show \1/ > 3 in at least one of the ^ values. Examples of 
the intensity profiles of high and low significance loops 
are shown in Fig. A.l. 

— Color index: AI 60 / AI 100 color indices were calcu- 
lated for our loops from the radially avergaged 60 and 
100 /im surface brightness profile. The color index is 
defined as the slope of the Iqo versus Iioo scatter plot 
using the data points of the surface brightness profiles 
in the positions of the loop wall only (a- m < a < a ou t). 
The uncertainties of AIqq / A/ioo given in the catalogue 
are the formal errors of the slope fitting. 

— Associated objects: We attributed associated ob- 
jects to our identified loops with the condition that 
an associated object has to be placed in the wall or 
in the interior of the loop. We considered the follow- 
ing type of possible associated objects (references are 
indicated): 

• dark clouds (Dutra & Bica l2"UfJ2"|) : 

• supernova remnants (Green 1994); 

• OB-associations (Lang 1992); 

• pulsars (Taylor et al. 1993); 

• HII regions (Sharpless 1959); 

• IRAS point source with molecular core FIR colours; 



• IRAS point source with T Tau star-like FIR 
colours. 

Selection criteria for molecular cores and T Tauri stars 
were selected from the IRAS Point Source Catalogue 
(Joint IRAS Science Working Group 1988) according 
to the following criteria: 

• point sources associated with galaxies were ex- 
cluded; 

• photometric qualities are 2 or better at 12, 25 and 
60 /im; 

• o molecular cores: 

0.4 < log 10 (!£) < 1.0 & 0.4 < log 10 (fjfi) < 1.3 
o T Tauri stars: 

0.0 < log 10 (fg) < 0.5 & -0.2 < log 10 (Sa) < 
0.4 

following the definitions by Emerson (1998). F 12 , F 2 § 
and Fqq are the 12, 25 and 60 /im uncorrected IRAS 
fluxes, respectively. 

List of Objects 

Ceoheus Bubble! ' on page 2 
Local Bubble! ' on page 4 
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Fig. 4. Galactic latitude distribution of loops in the 
outer (B + D, top panel) and inner (A + C, bottom panel) 
Galaxy. See Table 1. The binsize is 10°. 
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Fig. 3. Galactic distribution of loops with known dis- 
tances around the Solar System (located at the panel cen- 
tres). Top panel: location of infrared loops projected to 
the Galactic plane, up to a distance of <i = 5kpc. Bottom 
panel: magnification of the central region (d< 1 kpc). The 
distance from the Galactic midplane is indicated by col- 
ors as given by the color bars. The approximate po- 
sition of Spiral Arms (grey areas) were adapted from 
Gcorgclin & Georgelin (1976) The immediate outer sur- 
rounding of the Solar System was defined between I = 70°- 
250°, while the inner Galaxy space lies on the other side 
of the I = 70°-250° axis (dashed red line, see also Table 1). 
The majority of the loops are probably located in the 
Local Arm, as explained in Sect. 4.4). 
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Fig. 5. Average effective loop diameter vs. Galactic 
longitude in the inner (A + C) and outer (B + D) 
Galaxy; solid and dashed lines correspond to |6| < 30°and 
30° < \b\ < 60°, respectively. The binsize is 20°. 
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Fig. 6. (a) Distribution of estimated distances for 73 Fig. 7. Correlation of distance and effective size in the 
loops. Dash-dotted line represents the distribution of dis- "inner" (a) and in the "outer" Galaxy (b), based on the 
tances projected to the Galactic plane, (b) Distribution of sample of 73 loops with known distances, 
effective diameters in the same sample. The widths of the 
bins are 300 pc and 20 pc for (a, b), respectively. 
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Fig. 8. Location of known distance loops in regions A and 
B in the close vicinity of the Solar System. The size of 
the filled circles indicates the real effective size of the 
loops. Dashed and solid blue curves correspond to the 
W(D2)=20mA and W(D2) = 50mA co ntours of equiv- 
alent width of the Nal doublet line ( see ILallement et al.l 
2003). 
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Fig. 9. Distribution of effective diameters for known- 
distance loops on a logarithmic scale in the "inner" (a) 
and "outer" (b) Galaxy. The relation of log2i? ft-logA^ was 
fitted for 1.4 < log(2i? off ) < 2.6. The points used for the fit 
are marked by asterisks (see Appendix B). 
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Fig. A.l. Example of a (a) high significance loop (G077- 
77, ^noo = 38.5) and a (b) low significance loop (G019- 
19, ^noo = 3.7). Besides low surface brightness enhance- 
ments, bright spot-like objects and point sources in the 
neighbourhood of the loop may also result in poor sig- 
nificance values as was demonstrated in KMT04. Figures 
show the intensity profile before (top) and after (bottom) 
background removal (grey solid line) with the main de- 
rived parameters (J ex , <7 ex , a in , a out ). 



